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ABSTRACT 

As the largest gravitationally bound objects in the universe, clusters of galaxies may contain a fair sample 
of the baryonic mass fraction of the universe. Since the gas mass fraction from the hot ICM is believed to be 
constant in time, the value of the cosmological deceleration parameter qo can be determined, in principle, by 
comparing the calculated gas mass fraction in nearby and distant clusters (Pen 1997). To test the potential of this 
method, we compare the gas fractions derived for a sample of luminous (Lx > 10 45 erg s -1 ), nearby clusters with 
those calculated for eight luminous, distant (0.3 < z < 0.6) clusters using ASCA and ROSAT observations. For 
consistency, we evaluate the gas mass fraction at a fixed physical radius of 1 h$Q Mpc (assuming qo=0.0). We find 
a best fit value of go = 0.07 with -0.47 < qo < 0.67 at 95% confidence. This analysis includes both measurement 
errors and an intrinsic 25% scatter in the gas fractions due to the effects of cooling flows and mergers. We also 
determine the gas fraction using the method of Evrard, Metzler, & Navarro (1997) to find the total mass within 
rsoo, the radius where the mean overdensity of matter is 500 times the critical density. In simulations, this method 
reduces the scatter in the determination of gravitational mass without biasing the mean. We find that it also 
reduces the scatter in actual observations for nearby clusters, but not as much as simulations suggest. Using this 
method, the best fit value is qo = 0.04 with -0.50 < qo < 0.64. The excellent agreement between these two methods 
suggests that this may be a useful technique for determining qo- The constraints on qo should improve as more 
distant clusters are studied and precise temperature profiles are measured to large radii. 

Subject headings: galaxies: clusters of — galaxies: evolution — intergalactic medium — X-rays: galaxies — 
cosmology 



1. INTRODUCTION 

One of the major goals of cosmology has been to measure ac- 
curately the cosmological deceleration parameter qo- Determin- 
ing qo from the count magnitude relation for galaxies is notori- 
ously difficult due to evolutionary effects (e.g., Sandage 1988). 
However, recent studies of galaxy counts from redshift surveys 
provide reasonable constraints on qo . Schuecker et al. (1998) 
find a best fit of qo = 0.1 with an upper limit of qo < 0.75 at 
95% confidence assuming A = 0. Studies of high-redshift Type 
la supernovae (Perlmutter et al. 1997, Garnavich et al. 1998, 
Riess et al. 1998) suggest that qo may be negative due to a non- 
zero cosmological constant. Finally, a study of gravitational 
lensing statistics suggests that qo > -2 (Park & Gott 1998). 
These constraints do not, however, allow us to rule out any of 
the most popular cosmological models (e.g., closed universes 
with qo > 0.5 or f2 + A = 1 universes with qo ~ -0.5). It has 
been proposed that a constant Lx~T relation would allow cos- 
mological parameters to be measured (Henry & Tucker 1979). 
Another interesting method for improving or reinforcing these 
constraints is to analyze the gas mass fractions of rich galaxy 
clusters. Since this gas mass fraction is likely to be constant 
and its calculated value depends on cosmological parameters, 
comparisons of gas mass fractions of nearby and distant clus- 
ters can constrain qo (Sasaki 1996, Pen 1997). In this paper, 
we test the possibility of using the constancy of the gas mass 
fraction to constrain qo- We use the available data with suitable 
analyses to evaluate plausible uncertainties to determine if the 
method is sufficiently sensitive to allow accurate measurements 
of qo using future AXAF and XMM observations of large clus- 
ter samples. 

The approach we follow to constrain qo requires that the 



mean cluster gas mass fraction be a universal constant, indepen- 
dent of redshift. Rich clusters have masses ~ 1O 15 M , which 
suggests that the matter they contain ought to be a fair sam- 
ple of the distribution of baryons and dark matter in the uni- 
verse. Detailed analyses for hot, massive clusters (White et al. 
1993, Evrard, Metzler, & Navarro 1996; hereafter EMN) show 
that there is no known mechanism to significantly concentrate 
baryons (both galaxies and X-ray gas) in the center of the clus- 
ter or to expel them from the potential. This is supported by nu- 
merical simulations (e.g., White et al. 1993, Lubin et al. 1996), 
which show that the X-ray gas mass fraction f g = M gas /M, , 
averaged within ~ 1 Mpc should be constant, independent of 
redshift or evolution (Danos & Pen 1998, Frenk et al. 1998). 
In a flat universe, there is a degeneracy between redshift and 
a%, or equivalently, epoch and mass. Evrard (1997) has shown 
that the gas mass fraction is independent of total cluster mass 
to an accuracy of 1.5%. This is equivalent to an independence 
of the gas mass fraction in redshift to this accuracy. Alterna- 
tively, in an extremely open universe, clusters formation ceases 
early in the evolution of the universe. Afterwards, there is lit- 
tle evolution imposed on clusters by their environments. Most 
simulations (Thomas & Couchmann 1992, Katz, Hernquist, & 
Weinberg 1992, Cen & Ostriker 1993, Kang et al. 1994, Lu- 
bin et al. 1996, Pen 1998) show the gas to be slightly more 
extended than the dark matter within ~ 1 Mpc due to shock 
heating and other transfers of energy. This raises a possible 
concern, namely that the gas fraction may depend on radius. At 
the virial radius, inflow and outflow of mass will affect hot gas 
and dark matter identically, so that the gas mass fraction should 
remain constant at the virial radius. However, since the X-ray 
data do not extend to the virial radius, we determine the gas 
mass fraction at two radii - a fixed physical radius of 1 li^ Mpc 
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and a fixed overdensity radius rsoo (described below). If there 
are any large systematic variations of gas mass fraction with 
radius, these two methods should yield significantly different 
values of qo. As described below, they agree remarkably well. 

Evolutionary effects due to stripping of gas from galaxies 
and ejection of mass by supernovae should be relatively small. 
In nearby late-type galaxies, neutral gas contains <~ 10% of the 
luminous stellar mass. Since mass-to-light ratios of clusters 
are typically 300 h^, the gas associated with the galaxies in a 
rich cluster comprises ~3x lO -4 ^ of the total cluster mass. 

-3/2 

Since X-ray gas mass fractions are typically 0. 18/i 50 , this lm- 

1/2 

plies that the ratio of galactic gas to X-ray gas is <~ 0.002/z 50 . 
For early-type galaxies, Faber & Gallagher (1976) find that the 
amount of gas ejected into the ISM is 0.015 Moyr^lO 9 /,©)" 1 . 
For 1000 bright galaxies (L ~ lO n L0), this amounts to ~ 
2 x 10 12 M Q over a Hubble time (10 10 yrs) for the entire clus- 
ter. Since the mass in hot gas within the central 1 Mpc 
is ~ 1O 14 M for the clusters in the current study, this effect 
amounts to less than 2% of the hot gas mass in the core of the 
cluster. Thus, the effect of gas stripping from galaxies, both 
for late-type and early-type dominated populations, is almost 
certainly less than one percent of the mass in X-ray gas. 

Similarly, the ratio of luminous stellar mass to X-ray gas 
mass in clusters is 0. 06/159. Since only stars larger than about 
6Mq become supernovae, and since these massive stars com- 
prise a small fraction of luminous stellar mass, any injection of 
gas into the ICM from supernovae will be an effect of less than 
one percent. Thus, the actual gas mass fraction should vary by 
no more than two percent due to the effects of gas stripping of 
galaxies and ejection of gas from galaxies into the ICM by su- 
pernovae. Since the true gas mass fraction is assumed to remain 
constant while the calculated gas mass fraction depends on cos- 
mological parameters, any apparent change in the cluster gas 
mass fraction with redshift can be used to directly measure the 
evolution of the angular diameter distance (Sasaki 1996, Pen 
1997). As studied here, the cosmological dependence of f g is 

3 /2 

roughly f g oc d L , where di is the luminosity distance. The 
cluster gas mass fraction method provides a direct measurement 
of q {) (Pen 1997). The difference between a flat universe and an 
empty universe at a redshift of ~ 0.5 is 15%, significantly larger 
than the systematic effects described above. 

There have been some attempts to determine cosmological 
parameters using gas fractions from the literature (Sasaki 1996, 
Pen 1997, Danos & Pen 1998, Cooray 1998), but the gas frac- 
tions for distant clusters are found by different authors with 
unknown systematic differences in the measured fractions de- 
pending on the methods of data analysis. These studies also 
rely on small numbers of high redshift clusters. Since simula- 
tions predict an intrinsic scatter in the observed gas fraction of 
25% (Pen 1998, Danos & Pen 1998), a large sample is needed to 
avoid excessive noise due to small sample size. We present here 
the first "large" sample of cluster gas fractions at high redshifts 
determined using an internally consistent method of analysis. 
This study is also the first to select nearby clusters on the basis 
of X-ray luminosity, which is known to be strongly correlated 
with other cluster properties, most notably the gas temperature 
(e.g., David et al. 1993). Our sample consists of eight clusters 
at z ;> 0.3, which we compare to a sample of nearby clusters 
with comparable luminosities (L x > 10 45 erg s" 1 ). These are 
all of the distant (z ;> 0.3) clusters available in both the ASCA 
public archive and the ROSAT public archive which were ob- 



served with the PSPC and can be well-described by the standard 
hydrostatic-isothermal (3 model. Since the cosmological effects 
depend only on the apparent change of gas fraction, they are 
independent of the numerical value of the gas fraction. Thus, 
systematic uncertainties in the numerical value of the gas frac- 
tion (e.g., from uncertainties in Ho) are unimportant provided 
they are consistently applied to all clusters in the sample. 

To determine the total gas mass and total gravitational mass 
in a cluster, we must determine its temperature, temperature 
distribution, and gas density distribution. ASCA has a large 
energy bandwidth (0.5-10.0 keV), making it especially useful 
for measuring cluster temperatures, which are typically in the 
range 1-10 keV. However, ASCA has insufficient spatial reso- 
lution to accurately determine the gas distribution, so we must 
use ROSAT, with its higher angular resolution but lower energy 
bandwidth, to derive the gas distributions. As described below, 
we obtain -0.47 < q a < 0.74 and -0.50 < q < 0.64 at 95% 
confidence using two similar techniques for calculating the gas 
fractions. We describe the techniques for calculating gas frac- 
tions in Section 2 and the observations in Section 3. We calcu- 
late the cosmological dependence of the gas fraction and go in 
Section 4 and discuss our results in Section 5. 

2. METHODS OF CALCULATING GAS MASS FRACTIONS 

Since the dynamical time of the gas is much shorter than 
the age of the cluster, it is reasonable to assume that the gas 
is in hydrostatic equilibrium near the center of the cluster. In- 
deed, Jones & Forman (1984, 1998, hereafter JF98) found that 
a spherically symmetric hydrostatic-isothermal (3 model (Cava- 
liere & Fusco-Femiano 1976) is a surprisingly accurate descrip- 
tion of cluster profiles. This model provides good fits to the 
present data, so we use it in our analysis. The surface bright- 
ness profile, S(r), of a cluster can be accurately described as a 
function of radius r by a (3 model 



E(r) = £ 
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where a is the core radius of the cluster, (3 = fim p a 2 /3kT gas is 
the ratio of energy per unit mass in galaxies to the energy per 
unit mass in gas where a is the velocity dispersion, /1 is the 
mean molecular weight, and T gas (hereafter T or Tx) is the tem- 
perature of the hot gas (Forman & Jones 1982). j3 has been 
found to have a numerical value of 0.6 ±0.1 for 60 of 85 clus- 
ters (Jones & Forman 1998). This is different than the value 
of (3 ~ 0.8 found in some simulations (Evrard 1990), but this 
apparent discrepancy was resolved by Navarro, Frenk, & White 
(1995) as an effect of observational constraints. They further 
argue that the cause of the other "beta-discrepancy", the dif- 
ference between the value of (3 spec ~ 0.9 measured by spec- 
troscopy of galaxies (e.g., Bahcall & Lubin 1994) and the above 
value is due to galaxies not being distributed like an isother- 
mal sphere, as is assumed by the model. Recent ROSAT ob- 
servations of nearby clusters suggest that the surface brightness 
distribution of clusters is indeed accurately described by a (3 
model even out to large radii (Buote & Canizares 1996). For an 
isothermal gas, this yields a gas density profile described by 



^W = p (l + (^ z ) 



2 r 3/3/2 



(2) 



where r here (and for the rest of the paper) refers to the three- 
dimensional radius. The central density can be found from the 
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luminosity L(R) within a projected radius R by the relation 



L(R) = 



27rn e n//Aofl 3 
(1-3/3) 



l+s 2 + 







-3/3+1 



d+s 2 r 3l3+l 



ds 



(3) 

where Ao is the radiative cooling coefficient at the gas temper- 
ature, L(R) is the luminosity inside a radius R, s is the distance 
along the line of sight in units of core radius, and n e and 
are the electron and hydrogen number densities of the ICM at 
the center of the cluster (e.g., David et al. 1990). The total gas 
mass within a radius r is then just the integral of the density dis- 
tribution (Equation^). For ROSAT observations of hot clusters, 
David et al. (1990) find that the gas mass is fairly insensitive to 
a temperature gradient T(r). 

With the assumptions of hydrostatic equilibrium and spheri- 
cal symmetry, the gravitational mass inside a radius r is given 
by 

M tot «r) = -tt^ + -rr—jr (4) 

Assuming a uniform temperature distribution, the second term 
on the right hand side vanishes. We then only need to determine 
the gas temperature and the density distribution of the gas to 
calculate the gravitational mass. Combining equations^ and^j, 
we find that the mass is related to f3 and a by 



M tot {< r) : 



3kT(3r 3 



fim p Ga 2 (l + (^ 2 ) 



1.13 x io I4 /3r keV 



© 
(5) 

where M, ol (< r) is the total gravitational mass within a radius r 
and the numerical approximation is valid for T^y in keV and 
r and a in Mpc. The gas fraction f g (r) is then simply M gas /M lo , 
evaluated at a given radius. 

The assumptions of hydrostatic equilibrium, spherical sym- 
metry, and isothermality are not perfect. We know that the latter 
two assumptions are violated for some nearby clusters (Marke- 
vitch et al. 1998, Henry & Briel 1995), and clusters undergoing 
major mergers violate hydrostatic equilibrium. In fact, Marke- 
vitch et al. suggest that all clusters have temperature structure. 
The effect of these violations of spherical symmetry, isother- 
mality, and hydrostatic equilibrium is to introduce an additional 
scatter in the observed gas fractions beyond simple measure- 
ment errors. This scatter has been found by simulations (Danos 
& Pen 1998) to be about 25%. We include this intrinsic scatter 
in our determinations of qo in Section 4. 

Markevitch & Vikhlinin (1997) studied the variation of gas 
fraction with radius in detail for A2256. If A2256 is assumed to 
be isothermal, the gas fraction varies by only about 10% with 
radius. However, by including the temperature profile of the 
cluster, they found that the isothermal model overestimated f g 
at small radii and underestimated it at large radii, with the esti- 
mate being very close to accurate at a radius of 1 .2 Mpc. We 
assume an isothermal model with a radius of 1 /zjg Mpc, which 
overestimates the gas fraction by about 15% depending on the 
model used for the dark matter distribution. Similar results are 
obtained for a study of temperature profiles of 26 nearby clus- 
ters (Markevitch et al. 1998). More importantly, the shapes of 
the temperature profiles for these clusters was found to be ex- 
tremely similar. This similarity of cluster temperature profiles 
suggests that cluster-to-cluster variations in temperature profile 
will have only a small effect on the apparent change in gas frac- 
tion with redshift. Thus, since we lack detailed temperature 
profiles for distant clusters, one of our approaches is to assume 
isothermality and evaluate the gas fraction at 1 Mpc. 



The measurement of cluster gas fractions was explored in de- 
tail by EMN (1996) based on cosmological gas dynamic sim- 
ulations for several cosmological models. They find that clus- 
ter mass estimates based on the above method are surprisingly 
unbiased, yielding an average estimated-to-true mass ratio of 
1.02 with a standard deviation of 14%-29%. EMN also de- 
velop a method to reduce the scatter in the estimate of cluster 
mass based solely on the emission- weighted gas temperature Tx 
within a radius where the mean density is 500 times the critical 
density. This radius, denoted by r^c, is found to vary with Tx 
as 

Tx \i/2. 



r SO o = (2.48 ±0.17) 



lOkeV 



/z 5 oMpc. 



The mass within rsoo can be estimated as 



M 50 o(7i) = 2.22x 10 1 



lOkeV 



3/2 



Mr. 



(6) 



(7) 



and is found to have an average estimated-to-true mass ratio 
of 1.00 with a standard deviation of only 8%-15%. In Sec- 
tion 4, we use this method combined with the above method of 
calculating the gas mass in an attempt to reduce the scatter in 
the gas fraction and hence improve the constraints on qo- For 
A2256 (Markevitch & Vikhlinin 1997), the gas fraction at r 50 o 
is underestimated by 20-40% depending on the dark matter dis- 
tribution model. Again, this should mainly affect the absolute 
value of the gas fraction, not the apparent change in gas fraction 
as a function of redshift and thus should not significantly affect 
our results beyond contributing to the intrinsic scatter discussed 
above. 

Our use of two methods to determine the gas fraction can 
identify possible systematic effects due to changes in gas frac- 
tion with radius from observations. By comparing determina- 
tions of qo from cluster gas mass fractions evaluated both at a 
fixed physical radius of 1 Mpc and at a fixed overdensity 
radius (rsoo), we should be able to identify large systematic ef- 
fects caused by our inability to evaluate the gas mass fraction at 
the virial radius. 

3. OBSERVATIONS 

ASCA's broad energy band (0.5-10.0keV) is particularly use- 
ful for determining hot (> 5 keV) cluster temperatures. We 
therefore used ASCA data to measure cluster temperatures and 
luminosities, while we used the higher angular resolution of 
ROSAT to determine the parameters of the gas distribution. Be- 
cause of the poor angular resolution of ASCA, we determine the 
emission weighted average temperature within 6', which corre- 
sponds to ~ 2/ig, 1 ) Mpc at z=0.3. 

The temperatures and iron abundances of our sample are 
shown in Table 1. Details of the fitting techniques can be found 
in Rines et al. (1998). The temperatures we measure are consis- 
tent with those found by Mushotzky & Scharf (1997) for clus- 
ters common to both samples. The flux is measured between 
0.01 and 100.0 keV from the best-fit model, so it is effectively 
the bolometric flux. The bolometric luminosity is calculated 
from the bolometric flux of the best-fit model to the ASCA data, 
using the luminosity distance di = cHq 1 [z+ j(1 - qo)z 2 ] an d as- 
suming qo = 0.0. 

We used ROSAT PSPC observations for analyzing surface 
brightness profiles. Although the PSPC has less angular resolu- 
tion than the HRJ, its background is much lower, which is im- 
portant for extended faint objects such as clusters. To account 
for the point spread function of the PSPC, which depends on 
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Table 1 



Temperatures and iron abundances of several high redshift clusters from ASCA data. 



Cluster 


z 


kT 


90% 


Abund. 


90% 


Flux(10" 12 ) 


L£(boiy 


xi 


Counts" 






(keV) 


(keV) 






ergs cm" 2 s _1 


ergs s" 1 


^ ~ 1550 


(thousands) 


Zw3146 


0.291 


7.1 


6.7-7.5 


0.23 


0.18-0.28 


16.8 


7.99 


1.006 


37.6 


A1576 


0.302 


8.7 


7.5-10.1 


0.16 


0.04-0.28 


6.2 


3.22 


0.946 


11.7 


A2744 


0.308 


9.5 


8.8-10.6 


0.20 


0.11-0.29 


12.3 


6.79 


0.983 


21.8 


MS2137.3-2353 


0.313 


5.3 


4.8-5.8 


0.40 


0.29-0.51 


4.9 


2.76 


0.948 


7.2 


MS1358.4+6245 


0.328 


5.3 


4.7-6.0 


0.20 


0.09-0.32 


3.2 


2.00 


1.001 


7.5 


A959 


0.353 


5.3 


4.4-6.4 


0.05 


<0.28 


2.7 


2.00 


0.977 


4.6 


MS045 1.6-0305 


0.541 


9.0 


7.5-10.8 


0.19 


0.07-0.33 


3.2 


6.46 


0.975 


5.9 


CL0016+16 


0.555 


8.6 


7.3-10.2 


0.11 


<0.22 


4.0 


8.56 


0.976 


6.9 



"Luminosities are calculated assuming c/q = 0.0. 
'Total number of counts for all four detectors. 



the angle from the center of the instrument, we calculated the 
point spread function for each cluster depending on its location 
in the field of view (Markevitch & Vikhlinin 1997). The surface 
brightness distribution was fit to a (3 model convolved with the 
PSF with /3, £o, a, and a constant background as free param- 
eters using the method of maximum-likelihood. We exclude 
from our sample all clusters which cannot be fit to a physically 
reasonable (3 model (e.g., A851 and A1758, which are known 
to have significant substructure - see Schindler & Wambsganss 
1996, Rines et al. 1998). The fit for A1576 is shown in Figure 
1 , and the resulting parameters for all clusters are given in Ta- 
ble 2, along with the resulting gas masses, gravitational masses, 
and gas fractions as calculated assuming qo = 0.0. The param- 
eters for CL0016+16 agree with those found by Neumann & 
Bohringer (1997) for HRI observations. Zw3146 is known to 
contain a large cooling flow (Allen et al. 1996), so the tempera- 
ture and hence the total mass is likely underestimated. The ob- 
served gas fraction for this cluster would then be overestimated. 
With improved angular resolution, we could exclude the cool- 
ing flow region and determine a more accurate temperature for 
the cluster. 

Table 2 shows the gas fractions both at 1 Mpc and at rsoo, 
and we compare these gas fractions in Figure 3. These two 
methods use equations q] and ^ respectively to estimate M, ol ; 
both use equations || and || to estimate the gas mass. The val- 
ues of fg(rsoo) are lower on average than f s (lh^Q Mpc), but the 
difference lies mostly within the fractional errors taken from 
the 90% confidence limits on the temperature. Evrard (1997) 
finds that in simulations, the gas fraction at r^oo is related to the 
gas fraction at r x , according to fg(r 500 ) = f g {r x ) x (r 500 /rx) n , 
where r\ = 0.13-0.17. The above discussion of the effect of tem- 
perature profiles on f g suggests that f g (rsoo) is understimated 
by assuming isothermality. This effect is in the sense of reduc- 
ing the difference between the simulations of Evrard and the 
present results. 

4. COMPARISON TO NEARBY CLUSTERS AND CONSTRAINTS ON q 

To determine the apparent change of the measured gas frac- 
tion, we must know the gas fraction at the present epoch. We 
do this by averaging the gas fractions found by JF98 for clusters 
observed with Einstein at z < 0. 1 with luminosities greater than 
10 45 erg s" 1 and accurate measurements of (3 and core radius a. 
Comparing Einstein to ASCA data introduces only a small error; 



A1576 



100 - 
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1000 
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FIG. 1. — Surface brightness profile for A1576 and the best-fit f3 model plus 
background. Errors shown are rms pixel-to-pixel variations. 



Henry (1997) found that the fluxes determined from ASCA for 
EMSS clusters agree very well with those found by Einstein. 
For ten distant clusters, the mean flux ratio is 0.998, with an 
rms scatter of 0.205. There is excellent agreement between the 
gas fractions calculated from the method of equations || and || 
and those in JF98, so the relative gas fractions of nearby versus 
distant clusters should be accurate. The luminosity selection 
ensures that the nearby clusters are similar to those we study at 
high redshift. For example, the well-known relation between 
X-ray luminosity and temperature (e.g., David et al. 1993, 
Mushotzky & Scharf 1997) indicates that clusters with com- 
parable luminosities also have comparable gas temperatures. It 
has been found that the physical properties of individual clus- 
ters such as iron abundance, luminosity-temperature relation, 
and X-ray temperature do not evolve significantly for redshifts 
up to z ~ 0.4 (Mushotzky & Loewenstein 1997, Mushotzky & 
Scharf 1997). Thus, both simulations and observations sug- 
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Table 2 



Gas distribution parameters from ROSAT surface brightness profiles. 



Cluster 


z 




a 

Qql Mpc) 


'"soo 
(hr 5 l Mpc) 


Gas Mass 
(1O 14 M ) 
(1 h£ Mpc) 


Total Mass 
(1O 14 M ) 

(1^50 M P C > 


f* 

(1^50 M P C ) 


Gas Mass 
(10 14 M Q ) 

foot)) 


Total Mass 
(10 14 M Q ) 
O500) 


fg 

O500) 


Zw3146 


0.291 


0.617 


0.079 


2.09 


1.25 


5.05 


0.247 


2.97 


13.3 


0.224 


A1576 


0.302 


0.566 


0.179 


2.31 


0.908 


5.39 


0.169 


3.04 


18.0 


0.169 


A2744 


0.308 


0.887 


0.654 


2.42 


1.49 


6.67 


0.224 


4.85 


20.6 


0.236 


MS2137 


0.313 


0.642 


0.054 


1.81 


0.642 


3.83 


0.167 


1.25 


8.57 


0.146 


MS1358 


0.328 


0.641 


0.112 


1.81 


0.700 


3.79 


0.185 


1.42 


8.57 


0.166 


A959 


0.353 


0.732 


0.563 


1.81 


0.858 


3.33 


0.258 


2.17 


8.57 


0.253 


MS0451 


0.541 


0.929 


0.335 


2.35 


1.32 


8.47 


0.156 


2.78 


18.9 


0.147 


CL0016 


0.555 


0.869 


0.427 


2.30 


1.67 


7.14 


0.234 


4.19 


17.7 


0.236 



NOTE. — Masses calculated assuming go=0.0. 



gest that the dependence of the gas fraction on cluster evolution 
should be small for our sample. Finally, it is worth noting that 
while the above luminosity selection ensures that the two sam- 
ples of clusters are similar, we are not required by the method 
to make such a selection. For instance, Evrard (1997) showed 
that the gas mass fraction should be a universal constant for 
rich clusters and even groups of galaxies. However, selecting 
our nearby sample by luminosity removes a possible system- 
atic effect (though we shall see below that it makes very little 
difference in the results). 

We find that within 1 Mpc, the gas fraction of the nearby 

JF98 subsample is f g = (0.186 ± 0.013)/z^ /2 , where the stan- 
dard deviation for a single cluster is 0.048. This is consistent 

with Evrard (1997), who finds f g (r 500 ) = (0.170 ± 0.008)^ /2 
for a sample of nearby clusters with no luminosity cutoff. The 
observed scatter in f g is 26%, in excellent agreement with sim- 
ulations (Danos & Pen 1998). The calculated gas fraction of 
a cluster has the cosmological dependence f g oc h 3 ^ 2 flb/flo- 
Since simulations suggest that Q^/Qo is constant for a given 
cluster, we should be able to directly measure the redshift evo- 
lution of the Hubble constant h(z). We can derive the cosmo- 
logical dependence of the cluster gas mass fraction by 

f g = M gas /M tol oc p gas r 3 /r oc d L dl /2 /d A = d L d 1 / 2 . (8) 

Since dud^ 2 = df/ 2 /(l+z), where d A is the angular diameter 
distance, the dependence of gas mass fraction on qo for any 
individual cluster is f g (qo)/f g (qo = 0) = [d L (q )/d L (q Q = 0)] 3 ' 2 
or 

f g (qo) = (( 1 - q + q Q Z + (qo ~ 1 ) V 2 ^+ l)/(?o(z + z7 2 )) ) ^ 

Xf g (q = 0) (9) 

(Pen 1997, Zombeck 1990) where f g (qo = 0) is the gas fraction 
calculated assuming qo = and f g (qo) is the gas fraction calcu- 
lated assuming a different value of go- Since we assume qo = 0.0 
throughout our analysis, we must derive the correct value of qo 
using the above equation. 

Figure 2 shows the gas fractions (calculated using qo = 
0.0) of all of the clusters studied versus redshift along with 
the predicted apparent evolution in gas fraction for qo = 0.5 
(Q = 1,A = 0), q = 0.0 (Q = 0,A = 0), and q = -1.0 
(Qo = 0, A = 1). Errors shown are 90% confidence limits from 



the temperature measurements, which are the dominant source 
of error. The solid squares are clusters noted in JF98 as "sin- 
gle" morphological type, whereas the open squares are clusters 
observed to contain significant substructure. The scatter in the 
"single" clusters is significantly smaller than for clusters which 
are known to violate the assumptions of hydrostatic equilibrium 
and isothermality. Because we currently lack adequate spatial 
resolution to make similar analyses of the distant clusters in 
our sample, we must include the full 25% scatter in f g to place 
meaningful constraints on qo. Since the nearby sample is not 
located at z = 0.0, we also apply a correction to the nearby gas 
fraction assuming z = 0.06 (the mean value) for all of the nearby 
clusters. We define x 2 = £[/g(<7o,Zi) _ /g(<7o)] 2 / 7\ where of is 
the sum in quadrature of the 90% confidence limits in tempera- 
ture (a conservative estimate of the the l-a uncertainties in the 
observations), the 7% uncertainty in the present gas mass frac- 
tion, and the 25% scatter expected from simulations. With this 
definition, we find a best fit value of qo = 0.07 with an accept- 
able x 2 = 2.94 for 7 de grees of freedom. We calculate 95% 
confidence limits using A% 2 = 3.841 and assuming the errors 
are Gaussian. This gives a range of -0.47 < q Q < 0.67. If we 
assume the worst-case scenario of an actual evolution of 2% 
in gas fraction due to differences in gas stripping and injection 
from galaxies, the limits relax slightly to -0.51 < qo < 0.74. 

To apply the method of EMN to find the present value of 
fg( r 50o)< we must know the value of 7-500 an d the cluster tem- 
peratures, which we take from David et al. (1993). These pa- 
rameters and the resulting masses and gas fractions are shown 
in Table 3. We find an average gas fraction of fgir^oo) = 
0.176 ±0.011, where the standard deviation for a single cluster 
is 0.043 (the observed scatter is 24%), which is again consis- 
tent with the gas fraction of Evrard (1997) and the scatter in 
the simulations of Danos & Pen (1998). Our results support the 
hypothesis that the inferred gas mass fractions using the EMN 
estimator of total mass within rsoo yields gas mass fractions 
with slightly less scatter than those inferred using a constant 
radius of 1 Mpc, but the reduction in scatter is not as signif- 
icant as in simulations. Since the scatter for nearby clusters was 
not significantly reduced, we again assume an intrinsic scatter 
of 25% in the distant cluster gas fractions. Using the method of 
EMN, we obtain -0.50 < qo < 0.64 with 95% confidence (again 
assuming Gaussian errors). The best-fit is at qo = 0.04, with an 
acceptable x 2 = 3.90 for 7 degrees of freedom. If we assume 
that the scatter does actually decrease, the limits on qo narrow 
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Table 3 



Gas and gravitational masses and gas fractions for luminous nearby CLUSTERS (JF98) USING THE EMN 

METHOD. 



Cluster 


z 





a 


T x 


>"500 


Gas Mass 


Total Mass 


J s 








(hll Mnc) 


(keV) 


(hll Mnc) 

v 50 1TA FW 


(10 14 Mn) 


(10 14 Mn) 




A85 


0.0518 


0.62 


0.26 


6.2 


1.95 


2.03 


10.7 


0.189 


A401 


0.0748 


0.65 


0.26 


7.8 


2.19 


2.46 


15.3 


0.161 


A426 


0.0183 


0.58 


0.28 


5.9 


1.90 


2.85 


10.1 


0.283 


A478 


0.0881 


0.76 


0.30 


7.3 


2.12 


2.79 


13.9 


0.202 


A3266 


0.0594 


0.70 


0.55 


6.2 


1.95 


2.39 


10.8 


0.221 


A644 


0.0704 


0.70 


0.20 


6.9 


2.06 


1.53 


12.7 


0.120 


A754 


0.0528 


0.80 


0.58 


9.1 


2.37 


2.65 


19.3 


0.137 


A1650 


0.0845 


0.78 


0.29 


5.5 


1.84 


1.35 


9.06 


0.149 


A 1795 


0.0616 


0.73 


0.29 


5.6 


1.86 


1.72 


9.30 


0.185 


A2029 


0.0767 


0.69 


0.20 


7.8 


2.19 


2.50 


15.3 


0.164 


A2065 


0.0721 


0.64 


0.24 


8.4 


2.27 


1.84 


17.1 


0.108 


A2256 


0.0601 


0.73 


0.44 


7.4 


2.13 


2.31 


14.1 


0.163 


A2319 


0.0564 


0.69 


0.46 


9.9 


2.47 


4.01 


21.9 


0.183 


A3667 


0.0585 


0.51 


0.25 


6.5 


2.00 


2.40 


11.6 


0.206 



NOTE. — Masses calculated assuming qo=0-0- 




FIG. 2. — Observed gas fraction as a function of redshift, calculated assuming 
go = 0.0. From top to bottom, the lines show the apparent evolution of gas 
fraction predicted for q = 0.5 (f2 = 1,A = 0), qo = 0.0 (Cl = 0,A = 0), 
and q = -1.0 (flo = 0. A = 1). Errors shown are 90% confidence limits from 
temperatures and do not include the 25% scatter expected from simulations. 

slightly without changing the best fit value. Again assuming 
the worst-case scenario of 2% actual evolution in gas fraction 
due to galaxies, the limits become -0.53 < qo < 0.71. If we use 
Evrard's (1997) value of f g (rsoo) = 0. 170 for all nearby clusters 
(no luminosity selection), we obtain a best fit of qo = 0. 16 with 
95% confidence limits -0.47 < qo < 0.91, so the luminosity 
cutoff is not a source of significant systematic effects. Figure 
4 shows the variation of x 2 with qo for both methods. The ex- 
cellent agreement between the two methods suggests that there 
are no significant systematic differences. 



FIG. 3. — Gas fractions calculated from usual method (f g (\k^ Mpc)) and 
using EMN method to obtain M, ot (rsoo) and hence f g (rsoo)- Errors shown are 
90% confidence limits from temperature measurements. 

The low values of x 2 for the best fits above suggest that our 
error estimates are too conservative. A more robust way of 
calculating the confidence limits on qo is to use bootstrap re- 
sampling (Press et al. 1992). We drew 10 5 artificial data sets 
from the actual samples for each of our two methods. For the 
fixed physical radius method, 95% of the calculated values of 
qo lie between -0.18 and 0.66. For the EMN method, 95% of 
the calculated values of qo lie between -0.28 and 0.70. These 
calculations include none of the intrinsic 25% scatter from sim- 
ulations. If we again assume a maximum 2% actual evolution 
in the gas fraction, the limits change to -0.22 < qo < 0.74 and 
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FIG. 4. — x as a function of go- The solid line is the usual method of 
calculating gas fractions ( 1 h ^ Mpc) and the dashed line is using the method of 

EMN (/soo). The horizontal lines show 95% confidence levels (A% 2 = 3.841). 
The best fit for the usual method is go = 0.07, with 95% confidence limits of 
-0.47 < q < 0.67. The EMN method yields a best fit of 0.04 and limits of 
-0.50 < q < 0.64. 



-0.32 < qo < 0.78. The agreement with the above results is 
encouraging, as are the narrower confidence limits. 

5. CONCLUSIONS 

We have used the apparent change of the cluster gas mass 
fraction f„ = M gas /M, , with redshift together with the assump- 
tion of constant true gas fraction to derive constraints on the 
cosmological deceleration parameter qo- Our sample is the first 
"large" sample of high redshift clusters with a common method 
of analysis, which is important for an accurate measurement 
of the apparent change in f g with redshift. We compare high 
redshift clusters to a similar sample of nearby clusters, namely 
those with Lx > 10 45 erg s" 1 , though our results are not strongly 
affected by this selection. We include the intrinsic 25% scat- 
ter in f g predicted by simulations and verified by observations. 
We obtain -0.47 < qo < 0.67 at 95% confidence using the gas 
fractions at a constant physical radius of 1 Mpc and -0.50 
< qo < 0.64 with 95% confidence for gas fractions evaluated 
at a constant overdensity radius. The excellent agreement be- 
tween these two results suggest that this method to calculate qo 
is quite promising. Bootstrap resampling yields narrower confi- 
dence limits, -0.18 < qo < 0.66 and -0.28 < qo < 0.70 respec- 
tively for the two methods. Our results are in good agreement 
with the upper limit qo <0.75 from galaxy counts in redshift 
surveys (Schuecker et al. 1998) and the lower limit qo > -2 
from studies of gravitational lensing (Park & Gott 1997). 

We find that the standard closed, flat, and open models for 
the universe are all within the formal 95% confidence limits, 
as are some flat SI + A = 1 models with non-zero cosmological 
constant. We are able to rule out a universe with SI = and 
A = 1. Our results are consistent with observations of high- 
redshift Type la supernovae (Perlmutter et al. 1997, Garnavich 
et al. 1998, Riess et al. 1998) as well as with gravitational lens- 



ing statistics, which yield A < 0.65 (Kochanek 1996) at 95% 
confidence. 

Determining qo from distant cluster gas fractions has the sig- 
nificant advantage of being independent of Ho and the numer- 
ical value of f g . Systematic effects could be created by cal- 
ibration offsets, unknown selection effects, non-Gaussian er- 
rors, or cooling flows and mergers which are responsible for 
much of the intrinsic scatter in the gas fraction. We have in- 
cluded the expected intrinsic scatter in our calculations, and the 
results of Henry (1997) suggest that calibration offsets between 
fluxes from Einstein and ASCA are not very significant. The ac- 
tual cluster gas mass fractions may increase slightly with red- 
shift. This would be true if the gas in nearby clusters has been 
expelled from the core (e.g., by shock heating) systematically 
more than in distant clusters. Simulations suggest that this ef- 
fect is negligible in both S7 = 1 and SI < 1 universes. Alterna- 
tively, our results would be biased if the local value of Hq is 
significantly different from the global value, although "local" 
here refers to an average redshift of z=0.06. From Figure 2, we 
can see that the nearby gas mass fraction has much less scatter 
for "single" clusters (i.e., ones with no significant substructure). 
Thus, constraints on qo can be greatly improved not just by in- 
creasing the sample size but by improving the analysis of indi- 
vidual distant clusters so we can differentiate between clusters 
with and without substructure and properly measure tempera- 
ture distributions. In conclusion, although our current results 
do not improve significantly the determination of qo, we have 
shown that with the assumption of constant gas mass fractions, 
future detailed X-ray observations with AXAF and XMM can 
provide useful constraints on qo- 
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